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1 Introduction 



The properties of QCD at finite temperature have raised considerable interest in the 
hterature (see [|I| and references therein). At low temperatures it seems that color 
is confined and chiral symmetry is spontaneously broken. However, from asymptotic 
freedom, it is expected that at high temperatures, color will be liberated and chiral 
symmetry restored. It is a matter of intense debate whether there should be one or 
two phase transitions, at what temperatures they would occur, and what their nature 
would be. 

Within the standard wisdom the quark chiral condensate (0|gg|0) plays a central 
role in this problem, since it is assumed that chiral symmetry breaking is produced 
by a strong condensation of quark- ant iquark pairs [0]. However, in recent years, 
this hypothesis has been questioned, opening the possibility of small, even vanishing, 
(0|gg|0) scenarios |, H. (Note that we use (0|gg|0) for the condensate at T = 
and {qq) in general). 

The evolution of the quark condensate with the temperature has indeed been ad- 
dressed using several approaches. In general, the properties of {qq) can be derived 
from a somewhat idealized dilute pion gas, which is commonly described using an 
effective Lagrangian formalism [^], as we will do here, or by means of finite tempera- 
ture QCD sum rules 0. In general, all these and other approaches yield a rather 
consistent picture, although they usually have the large condensate assumption built 
in. 

In this work we want to know how the actual value of (0|gg|0), as well as the light 
quark masses, can modify the behavior of the chiral condensate, as for instance, with 
changes in the phase transition temperature. With that purpose, we will describe 
the pion gas by means of the virial expansion and using the interactions obtained 
from the Chiral Perturbation Theory (ChPT) formalism |0, although allowing for a 
wide range of (0|gg|0) values. Such a framework is usually referred to as Generalized 
Chiral Perturbation Theory (GChPT) |^ It should be noticed that we will be 
deahng with two effects: First, at T = 0, the size of (0|gg|0), which may be different 
from the standard large value. Second, the evolution at finite temperature which is 
also changed through the modifications in the meson interactions due to the different 
scheme of the explicit chiral symmetry breaking. Our purpose is to study what the 
interplay is of these two effects. 

The plan of the paper is as follows: In section one we describe briefly the ChPT 
and GChPT formalisms, focusing on the relation between the quark and meson masses 
with the quark condensate. The next section is devoted to the virial expansion for the 
pion gas, where we introduce the temperature dependence. In section three we make 
the study of the condensate dependence both on the temperature and the ratio of 
light quark masses, using the 0{j)^) amplitudes of GChPT. Next, in section four, we 
estimate the contributions from heavier states, and in the conclusions we summarize 
our results. 
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2 Standard and Generalized Chiral Perturbation 



Theory 

When considering just three massless quark flavors, the QCD Lagrangian exhibits an 
SU{3)l X SU{3)ii symmetry which, even neglecting particle masses, is not present 
in the physical spectrum. Instead, we observe an approximate SUl^r{3) symmetry, 
which means that the SU{3)l-r chiral group has to be spontaneously broken. Ac- 
cording to the Goldstone Theorem, there should be eight massless Goldstone Bosons 
(GB), which are identified with the pions, kaons and the eta. In a first approximation, 
these GB couple to the spontaneously broken currents with strength F ~ 90 MeV. 
These particles are so light compared with the typical hadronic scales, that they will 
dominate the hadronic dynamics at low energies or temperatures. 

In order to describe the hadronic dynamics at low energies we can therefore use 
these fields to build an effective Lagrangian, made of the most general terms that 
respect the above symmetry breaking pattern. As we are interested in the low energy 
regime, the terms are organized according to their number of derivatives. It can be 
seen, by counting the powers of momenta of different diagrams, that it is possible to 
renormalize any calculation and obtain finite results order by order in the expansion 
T0| . We could also couple gauge fields, scalar and pseudoscalar sources, etc., which 



would allow us to describe other processes. This whole approach is usually known as 
Chiral Perturbation Theory (ChPT) |^. 



Explicit Chiral Symmetry Breaking 

Up to the moment we have just considered the chiral limit. When quark masses are 
turned on, the GB become massive pseudo-GB and their masses can be obtained, 
generically, as 

M2~25om + O(mj) 
M|-~(m + ms)Bo + C(m^) 

M2~^(m + 2m,)5o + 0(mj) (1) 

where m = (m„ +771^)/ 2 (we will consider isospin as an exact symmetry) and Bq and 
other coefficients that may appear at higher orders are to be determined phenomeno- 
logically. Throughout this work, the first one will play a very relevant role, since it 
has a very physical meaning: In the chiral limit, and up to a normalization factor, it 
is nothing but the chiral condensate; namely 

{0\qq\0) = {0\uu + dd\0) " -2F^Bo (2) 

At this point two different approaches appear in the literature. The first one, still 
called ChPT 0, is to assume that the mass expansions in eq.(|I|) are dominated by the 
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Bq term. Its origin can be traced to the Gell-Mann-Okubo (GMO) and Gell-Mann- 
Oakes-Renner (GOR) formulae, which, within the effective formahsm, are obtained at 
first order by ehminating Bo in eqs.(|l|) and (Q). This large condensate scenario usually 
requires Bq ~ O(lGeV) and, apart from the GMO and GOR formulae, it is supported 
by several lattice calculations [|r^. Within this framework, the quark masses count 
as O(p^). The second approach, know as Generalized Chiral Perturbation Theory 
(GChPT) P, §], is nothing but considering the possibility that the (9(m^) terms 
could be of comparable size or even larger than the Bq term. As a consequence, 
both the quark masses and Bq count as 0{p). This approach is supported by some 



deviations from the Goldberger-Treiman relation in vriV, KA and KH [|T2| and some 
calculations using variationally improved perturbation theory or a relativistic many 
body approach |T3 . 



Those two alternatives are usually compared with the spontaneous magnetization 
M of spin systems: On the one hand, ferromagnets present an ordered ground state 
where the magnetization spontaneously acquires an M 7^ value. That would be 
analogous to the standard ChPT. On the other hand, in anti-ferromagnets the mag- 
netization remains at M = 0, which would be similar to the extreme case of GChPT 
where Bq = 0. Note that, despite their difference, in both systems the spins are ori- 
ented in one preferred spatial direction and therefore the 5*0(3) rotational symmetry 
is broken. 

Back to our subject, it should be noticed that both approaches have the same 
terms in the Lagrangian, although they are organized differently, and their relative 
size is also changed. Indeed, it is possible to reobtain standard ChPT as a special 
case of GChPT. 

At present, the experimental data does not exclude any of the two scenarios, 
although this question may be solved in a few years with an accurate measurement 
of TTTT scattering lengths from the decay of tt+tt^ atoms fl^ . 



Thus, since we are interested in high temperature differences with the standard 
scenario, throughout this paper we will use the GChPT formalism. As usual, the 
pseudo-GB fields are grouped in an SU(3) matrix as follows: 



U = exp{i^/F) ; $ = 72 



(3) 



And then, with the GChPT power counting described above, the Lagrangian 
is usually written as 

£(2) = i_ {tr(D^f/Z}^f/t) + 25otr(A<(f/t + U)) 

+Ao tT{MU^MU^ + MUMU) + ii{M{U + U^)Y 

+ ii{M{U - U^)f + 2Hoti{M'')] (4) 
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where A4 = diag{rh,rh,ms) is the quark mass matrix. In standard ChPT, only the 
two first terms are O(p^), whereas the rest is counted as O(p^). From the above 
Lagrangian we obtain the following meson masses 

Ml = 2rhBo + Atji^Aq + 4m(2m + m^)Zo 

M'^ = {m + ms)Bo + (m + nisfAo + 2(m + m^)(2m + ms)Z^ 

M^ = '^{m + 2m,)Bo + ^(m^ + 2ml) 

4 8 
+ -(m + 2m,)(2m + m,)Zl + -(m, - mfZ^ (5) 

Comparing with eq.(||), we have just added the Oirn?^) terms. In the standard for- 
malism, since only Bq is present, it can be eliminated and one recovers, at O(p^), 
the GMO and GOR relations. That is no longer possible in GChPT, although these 
relations will be recovered at higher orders. Of the three 0{m?^ parameters there 
are two, Z° and Zp, which violate the Zweig rule. They are expected to be small 
from large N^. arguments and is usual to neglect their contribution, and so we will do 
likewise in most of what follows. 

Note that, since the pion, kaon and eta mass values are known, then, changing the 
value of Bq is nothing but changing the values of the quark masses. As a matter of 
fact, the ChPT relations are frequently used in the literature to obtain ratios of light 
quark masses (for a recent update, see [l^ and references therein) or even to evaluate 



rh itself. However, most of these works have used the standard ChPT formalism and 
have the large condensate assumption built in, so that their results would change if it 
was removed. Nevertheless, there are determinations of irts — which do not rely 
on a large condensate value. For the sake of simplicity, and in order to facilitate the 
comparison with previous works Q], we will use the value — m„ = (184 ±32) MeV, 
given in [11^. That is, 

184±32_ 

m = MeV (6) 

r — 1 

(There are other similar analyses in ref.[|T^], whose results are all consistent with the 
previous relation.) As a consequence, the parameter that determines the relative size 
of the 0{mq) and (9(m^) terms is the quark mass ratio r = rris/rh, which ranges in 
the interval 

n^2^-l<r<2^-l^r, (7) 

The upper limit corresponds to the extreme case of a very large Bq condensate, 
whereas the second corresponds to i?o = 0. (Vacuum stability requires Bq,Aq,Zq > 
0). 

Of course, all these formulae are valid up to 0{p'^). For the moment, we have 
restricted ourselves to the OijP') case since it already displays the features of GChPT 
which are relevant for this work. In section four we will state our results including 
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Figure 1: The chiral condensate at zero temperature as a function of the hght quark 
mass rh, eq.(y). The continuous hne corresponds to the central value of Q. The 
uncertainty due to its error is covered by the shaded area. 



higher order corrections, although we will just present the GChPT formulae without 
such a detailed introduction. 

The chiral condensate at zero temperature 

Using the GChPT Lagrangian in eq.(P, the chiral condensate at 0{p^) is then given 
by 

(0|gg|0) = -2F^iB + m(Ao + Ho) + ...) (8) 

where B = BQ + 2{ms + 2m)ZQ . In practice Bq cannot be separated from B by looking 
at quark masses alone, but we have already commented that Zq is expected to be 
very small, so that B Bq. The parameter Hq is associated with the contact term 
of two scalar sources, which does not contain meson fields. However, it is needed as a 
short distance counterterm, and it indeed depends on the renormalization conventions, 
which introduce some small ambiguity (see |^ for a discussion). Nevertheless, using 
QCD sum rules with a simple model for the spectral function, and keeping F^M^ 
fixed at its physical value, it has been found [0 that the chiral condensate can be 
described by 

m{Q\qq\0) ^ -F^ [m^ - Am^n] (9) 

with = 4.7 ± 0.7. At 0{p^), the parameter is nothing but (Aq - Hq) /2. We will 
use the above equation to estimate the size of the quark condensate at T = 0. 

In Fig.l we show the dependence of the condensate with m, for Q = 4.7. Note 
that the plot starts aim = 1 MeV, which is approximately the standard ChPT value. 
There, (0|gg|0) ~ — (280)^MeV^, and it decreases smoothly as m gets larger, until 
it vanishes around m ~ 30 MeV. The shaded area between dashed lines covers the 
uncertainty in VL. 
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As a check of eq. we can see that it is consistent with previous estimates within 
the framework of standard ChPT where 

rn(0|gg|0) = -F^M^ 

and the hi and I3 parameters play a similar role as that of the GChPT pa- 
rameters, although in Standard ChPT they appear at 0{p'^). In that case, estimates 
based on a simple p resonance model and the large A^^^ limit yield c = 0.87 and 
c = 0.90 ±0.05, respectively 0. If we introduce in eq.@ a value in the range m from 
5 to 10 MeV, and take into account the fact that in standard ChPT Ft^/Fq ^ 1.057, 
we obtain c = 0.95 ± 0.4, which is a highly non-trivial check of eq.@. (Throughout 
this section we have neglected higher order logarithmic contributions that would yield 
corrections of the order of 1%) 



1 + 



M2 

TT 

327r2F2 



{4hi + h-r 



~F^,Ml- 



(10) 



3 The virial expansion and temperature effects. 

At low energies the free energy z is dominated by the contributions from the lightest 
particles. Therefore, we can use the Euclidean form of the above Lagrangian (denoted 
L{x)) within the standard finite temperature functional Euclidean formalism. Hence, 
in the thermodynamic limit, 

z = -Tlim^^oo^ / [dU] exp f - / d'^x L{x)\ (11) 
L-^ J \ JL»x[o,T] y 

where, as usual, the functional integration is over pion fields which are periodic in the 
Euclidean time, with period (3 = 1/T (see ref.Q). From the free energy we can derive 
any other thermodynamic property of our system, but let us first notice that since 
there is a spontaneously broken symmetry, even at T = there is some non-vanishing 
vacuum energy density eo- As a consequence, the pressure is defined only from the 
temperature dependent part of the free energy, P = eo — z. 

The quark condensate is now obtained as the derivative of the free energy with 
respect to the quark mass. That is 

m = 1^ = (0|gg|0) - ^ (12) 

where we have used that at T = the condensate is nothing but the vacuum expec- 
tation value (0|gg|0) = deo/drh. 

In this section we will just concentrate on how to obtain dP/ dm. For that purpose, 
one possibility is to calculate the free energy from the effective Lagrangian, as was 
done in within standard ChPT. That method follows the very same philosophy 
of the chiral expansion, but is rather lengthy. In this paper we will make use of 
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existing one loop calculations of elastic vrvr scattering, together with the relativistic 
virial expansion of a pion gas |T9|, ^ . 



Let us then consider a gas made only of pions. This approximation seems rea- 
sonable as long as the temperatures remain sufficiently below the kaon threshold [Q. 
Within the virial formalism, the pressure can be expanded as follows 

P = 3T(^) j:Bke-^''^' = 3^^j:BkiT)e. (13) 
\ ^ k=i ^ k=i 

The factors of three come from the fact that we are neglecting isospin breaking effects. 
Thus, there are effectively three different species of particles, labeled according to their 
third isospin component, that behave identically with respect to strong forces. We 
have also defined A = {2ti / M^^TY^'^ , which is the thermal pion wavelength. Note that 
the expansion parameter is the fugacity ^ = e~^^^ . In a non-relativistic framework, 
the expansion is usually performed using the definition ^ = e^'^, where /i would be the 
chemical potential. In contrast, in the relativistic case, there is a rest energy given 
by M^, whose contribution to eq.(|ll|) is equivalent to a chemical potential /i = — 
in a non-relativistic description. 

There is a closed expression for the virial coefficients of the free gas, which is 

In order to deal with the interacting gas, we will just consider two particle interactions, 
which can be justified as long as the density remains small. In it was shown that 
this is consistent with the three loop calculation in ChPT. In such case, it is enough 
to keep the two first terms of the virial expansion, whose coefficients will be given by 



B,{T) = Bf\T) (15) 



where Ki{x) is the modified Bessel function which behaves as y7r/2xe~^ when x — > 
oo. It is important to notice that the only dynamical information we need are the 
phase shifts 6ij. As an estimate of the applicability of this approach, it was shown in 
pO| that the second order virial expansion yields less than a 1% error when applied 
to the free gas up to T ~ 250 MeV. 

Let us finally remark that the high temperature behavior of the chiral condensate 
will be due to two different effects. First, the starting T = value, which may differ 
from the standard, large condensate, value. But, second, it also depends on how the 
mass dependence of the phase shifts has changed with respect to standard ChPT. 



8 



In the next section the phase shifts will we obtained using the existing GChPT 
calculations of the vrvr elastic scattering amplitudes [Q. In later sections we will 
include contributions from particles more massive than pions. 



4 The General Scenario 

Higher orders in GChPT 

Within the standard ChPT it was shown in that the O(p^) contributions accelerate 
the melting of the chiral condensate, lowering the critical temperature. Our aim now 
is to include the equivalent corrections within GChPT. Unfortunately, we have already 
seen that the Lagrangian has more terms that the standard C^'^\ That means 
that there are many more phenomenological parameters in the Lagrangian, which in 
many cases are not very well known. The situation gets even worse at higher orders. 
In general, the GChPT Lagrangian is built of terms like 0, |^ 

^^'^^ = B^'^^k^i), with £(,,0 ~ 0{p'm[) (16) 

k+l+n 

Indeed, we have already given 

£(2) 

in eq.(H) and we found that some of the constants 
are not very well determined. For the complete expression of the O(p^) Lagrangian 
we refer to ||^. For our purposes, there are several relevant modifications to our 
previous discussion: First, now there are three different decay constants F^, and 
Ffj. Second, neglecting Zweig rule violating parameters, the expressions for M^- and 
Mk in eq.d^) are now modified to 

=2mEo + 4m2Ao + ^5M^ 

3.Ml = {m + m,)5o + (m + nisfAo + (17) 

where 5Mf are higher order corrections and logarithmic terms , whose size is 5Mf < 
O.lMf ||2T[ (see below). As a consequence, the range of allowed r values is shifted 
upwards to 



,_^ FkMk . , ^JFkMk 
r, = 2 1 < r < 2 



,2 



With this modifications now ^ 8 and r2 can be as large as 39. 

Finally, there are also higher order corrections to the T = condensate itself, 
which contain chiral logarithms. That means that we cannot simply say that = 
{Aq + Hq)/2. Nevertheless we can still use the phenomenological parameter Q = 
4.7 ±0.7 in eq.(g). 
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The one loop tttt amplitude in GChPT 

Next, we need the tt+tt^ ti^tt^ scattering amplitude itself. Although it has been 
calculated in GChPT up to two loops |P, ^, for our purposes it will be more than 
enough to consider the one loop result, which reads: 



A{s, t, u) 



: + — 



3 ^ 



(19) 



+ - 2M-f + |l [it 



Ft 



2M. 



2\2 



[U 



2M\ 



2\2 



[s,t,u) 



where 



S, t, u) 



5 ^'2 ^ nf. 



+ 



12F4 
1 



+ 



12Fi 



2 / 4 ^ 2 

-aMl -f3[t- -Ml 



2-^2 nf. 4 2 



3"^^^-n^ 3- 



+ /?2(s-M)(t-4M2 



+ (3\s-t){u-AMl] 



J(t) 

J(n)(20) 



and J is the standard one-loop integral 0. 

In the literature, the values of the a, (3, Ai and A2 are fitted from experiment. 
However, in order to obtain the condensate dependence with the temperature, we 
need the derivative of the pressure with respect to M^r, and just a fitted value is not 
enough. Therefore, we also need to know the dependence of the parameters, and, 
if we want to study the effects of changing the light quark masses, we also need the 
dependence on r. 

Phenomenological parameters 

The actual expressions of the a, /?, parameters are rather complicated and involve 
many parameters from the GChPT Lagrangian, which frequently are not very well 
determined. In addition they contain chiral logarithms. It is therefore very convenient 
to expand a and (3 in powers of quark masses, namely 



a 



E 



a 



[n) 



/3 = E f^^'^- 



(21) 



n=0 ri=0 

Notice that, in GChPT, since quark masses are considered as 0{p), these expansions 
involve terms that count as odd powers of momenta. 

The above expansions have been worked out in 0], and they are the following: 



a(r) = l + 
(3{r) = l + 



— 1 
2 / 



' F2 



1 \F3 



1 + 18(2 -r)pi- 6rp2 + a^^\r) 



r-l\F3 



(22) 
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where in all the above equations we have neglected the Zweig rule violating parame- 
ters. 

Let us now try to estimate the size of the different terms in the a and /3 expansions. 
Let us then look back at the allowed values of r, eq.(|T8|). The relevant point for our 
discussion is that now, even with the lowest value r = 8, we obtain, using eq.(||), 
that ra < (26 ± 4.6) MeV. Therefore we can estimate that the terms coming from 
Bq and Aq should be C(l), those from C^^^ should be 0(10%) and those from £(2,2) 
and £(0,4) should at most reach the 1% level. Consequently, we will neglect the a^^^ 
and P^"^^ effects. The only parameters that remain undetermined are the pi^2, which 
contribute to a^^\ However, from a dimensional analysis [Q, their magnitude can be 
naively estimated as \pi\ ~ (0.4±0.2)/(r — 1)'^. Their dependence on the actual value 
of (needed for the numerical derivation) seems very weak. In our calculations we 
will take them first as zero and then we will include them in the uncertainty. 

Concerning Ai and A2, they come only from the terms in >C(4^o)) which do not 
contain explicit chiral symmetry breaking. They are given by 

A.^Af = 4(2UM + L3) - (log + i logM + ''^ 



487r2 [ fi^ 8 /^2 24 
1 



A.=Af = ,L'M - |,og-^ + -log:^ + (23) 

It can be seen that these parameters do not carry any r dependence. For definiteness, 
we will use for them the values obtained in 0]: 

Ai = (-5.3±2.5)10~^; A2 = (9.7 ± 1.0)10-^ (24) 

which are consistent with other determinations in the standard framework. 

The values of a and f3 depend on whether there is actually a large or small con- 
densate at T = 0, and we will use their r dependence to reproduce different scenarios. 
For illustrative purposes, let us recall that in the standard formalism both a and (3 
are slightly bigger than one and r ^ 26. In contrast, the low condensate alternative 
seems to prefer a ~ 2 and r ~ 10 

Phase shifts 

For the virial expansion we need the phase shifts of definite isospin and angular 



momentum channels. At lowest order, they are defined as follows (see |^2| for a 
discussion on this subject) 

tan 5,j(s) = (T(s)Re(t/j(s)), (25) 



where a{s) = yl — 4M^/s. The partial waves tj^j are obtained from the isospin 
amplitudes 

To(s, t, u)=3A{s, t, u) + A{t, s, u) + A{u, t, s), 
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Ti(s, t, u)=A{t, s, u) — A{u, t, s), 
T2{s,t,u)=A{t, s,u) + A{u,t, s), 



(26) 



by means of 

^iJ=^l d{cos9)Pj{cos9)Ti{s,t,u). (27) 

where Pj is the corresponding Legendre polynomial. In our calculations we have 
just used the lowest angular momentum for each isospin channel, namely (J, J) = 
(0,0), (1,1) and (2,0). For all means and purposes, they dominate the low energy 
pion interactions. 



The calculation of dP/drh 

We have then used the above phase shifts with the second order virial expansion. In 
order to obtain the condensate, eq.(|12|), we then need dP/drh, which can be obtained 
using 

dP dP dM^ dP SMk dP dF^ dP dFx , , 
_| _| _| 11 (28) 

dm OMt, dm OMk dm OFt, dm OFk dm 
Naively one just expects the first term, but let us remember that M^-, F^^ and F^ are 
m-dependent and they appear in the amplitude either directly or indirectly through 
a, /3, Ai and A2. That problem was carefully avoided in @] by using SU{2) standard 
ChPT and only using F in the free energy expansion. 

Of course, only appears in the fugacity, or in the free gas virial coefficients 
and thus we expect the three last terms in eq.(^) to be much smaller than the first. 
Indeed, within the range of r and T that we are interested in, we have found that the 
term due to the appearance of Mk in the amplitude is smaller than 1% and we have 
neglected it. In contrast, and Fk together generate contributions of the order of 
5%, and therefore they have been included in our calculations. 

The derivative of the pressure with respect to Mjr, -F^r and Fk have been performed 
numerically, with an increment of O.lMeV. For instance, the value of the pressure is 
first calculated with the real M^^ and then with — 5M^, including a change in the 
chiral parameters, following eqs. (p2D and (pSj). A similar procedure is followed for 
and Fk- 

In our calculations we have used 



1 + 2%^ 
— 1 



dm 2m 
OMk Ml r(2r* - r 



dm AuiMk — 1 

dF^ F^ (F, 



dm "'m[{r-l) + (F^/F^ - 1)] F^ 
dF^ r-1 dF^ 
dm 2 dm 



1 



(29) 
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Figure 2: Quark condensate versus temperature and r for = 4.7. 

which we have obtained from eqs. ([T7|) and from 0]. There are, of course, corrections, 
but their effects on the final resuhs are again less than 1%. 

Results 

As we have already commented, the virial expansion can be trusted only at low 
temperatures, mostly, due to the fact that above ~ 150 MeV the contributions from 
other more massive particles becomes relevant. These effects will be studied in the 
next section and we will see that they tend to lower the critical temperature, which 
is therefore more favorable for our approach. For the moment, if we give in our 
figures results for higher temperatures, they should be interpreted with great care as 
a qualitative behavior or, for instance, as a tendency towards symmetry restoration. 
Nevertheless, comparing between different figures could also illustrate what is the 
qualitative effect of a change in the parameters. 

Thus, in Fig. 2 we have shown the dependence of the chiral condensate with the 
temperature and for a light quark ratio in the range 8 < r < 26. For Q we have used 
the central value 4.7. Although the actual points at which (qq) = are just gross 
estimates, we can see that lowering r yields a systematic decrease in the chiral phase 
transition temperature. 

Indeed, Tc seems to be above 200 MeV for r > 20 going down to 130 MeV around 
r = 8. Note that for the latter temperature our approximations can become quite 
reliable. As we have already seen, smaller values of r are forbidden to ensure vacuum 
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Figure 3: a) Estimate of the errors in the {qq) evolution due to the uncertainty in f2. 
b) Uncertainties in the parameters that appear in the vrvr phase shifts, added hnearly. 

stabihty. 

The previous results have been obtained using the central values of all parameters. 
In Fig. 3a we show what happens if we take into account the uncertainty in Vt and 
thus, we plot the temperature dependence for the two extreme cases, r = 26 and 
r = 8. The former, which corresponds to the upper curve, is almost insensitive to 
this variation. It corresponds to the standard formalism, where the value of the chiral 
condensate is largely dominated by the 0{m) term and, consistently, changes in the 
other terms are almost negligible. On the other curve, which is associated with the 
lowest condensate scenario, the effect of this error is translated to a 10 to 15% change 
in T(7, at most. 

In Fig. 3b we show the uncertainties associated with all the parameters that 
appear in the scattering amplitude. In the shaded areas, we have taken into account 
all the effects of changing Ai, A2, pi and pi . In addition we have also let the 
pion and kaon masses vary between their values for the neutral or scalar particle. 
Note that in the case of the pion mass, such a change also affects the coefficients of 
the virial expansion and the fugacity. Finally, we have included the uncertainty in 
Fk/F^ = 1-22 ± 0.01 and we have let F^ change between 92.4 and 93.2 MeV which 
are two values currently cited in the literature. Both and do also appear in 
the expression of the T = chiral condensate. All in all, the overall uncertainty in 
Tc due to these parameters seems to be of the order of ±5 MeV at r = 26 and ±3 
MeV at r = 8. Since we have just simply added the different errors, we consider these 
numbers ClS db conservative estimate. 

5 Other massive particles 

In this section we will consider the effect of adding heavier particles to our pion gas. 
We will be following closely the approach of Gerber and Leutwyler ||^ with slight 
modifications to implement also the low condensate scenario. 
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The density of massive states should be exponentially suppressed by Boltzmann 
factors exp(— Mj/T), which means that their two body interactions will carry an 
exp [—{Mi + Mj) /T] factor. In addition, their interactions with pions are also sup- 
pressed by T'^/F'^, due to the chiral symmetry. Hence, we can treat those heavier 
particles in the free gas approximation. In such case, we have an additional contri- 
bution to the pressure, which is given by 



(30) 



where Qi is the state degeneracy of a state with mass Mj (that was the factor of 3 in 
eq(p!3D). Note that, since we will be dealing with temperatures much smaller than the 
first hadronic fermions, it makes sense just to use Bose statistics. The above formula 
is only meaningful at low temperatures, since as we increase the temperature, the 
mean distance between massive states shrinks and the dilute gas approximation is no 
longer valid. In ref. |^ it was estimated that the model is valid up to temperatures 
on the order of 150 MeV, although it "rapidly deteriorates" for higher temperatures. 

Back to the condensate, and in view of eq. ([T^) , the new contributions are of the 
form 

^ dMi dm ^ ' 



and therefore 



1 PiM. /-oo r)2 1 

Thus, we only have to estimate the value of dMi/drh. Naively, one would expect that 
the contribution m to a hadron mass would be roughly proportional to the number 
Ni of u and d quarks it contains. That estimate seemed quite appropriate in the 
standard framework [Q. We now have to check that it is also the case in GChPT. 

Let us then go back to eqs.(|^), since to get rough estimates it is enough to work 
at OijP'). As usual, we neglect the Zq and Zq parameters. Then, we obtain the 
following derivatives 

r(2r2 — r) — 1 





Ml r 


dm 


~AmMK 






dm 


GTTlMrf 



r2 - 1 
To — r 
— 1 



(33) 



We can reproduce the standard scenario with r = 26, which yields m ~ 7.4 ±1.3 
MeV using eq.(^). In such case, we find dM^/drh ^ 1.3 ± 0.2, which is in very good 
agreement with a rough estimate of 1. We also find dMr,/dm ~ 0.8 ± 0.1, again 
consistent with the naive estimate of 2/3. In any case, it seems that dMi/dm = Ni is 
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Figure 4: The evolution of the chiral condensate when we include corrections from 
a free gas of particles more massive than the pion. The shadowed regions cover the 
uncertainties in dMp/drh described in the text. These contributions always tend to 
lower the critical temperature. 

a small underestimate of the actual values of the standard scenario, as it was already 
pointed out in 0. Thus, in that work they considered that the range from Ni to 2Ni 
was a "fair representation" of the uncertainty in dMi/dm. 

However, if we set r = 8, which corresponds to the lowest allowed T = conden- 
sate, we find dMK/drh ^ 2.0 ± 0.4 and dMr^/drh ~ 0.36 ± 0.06. Again, the order of 
magnitude is correct, although within a factor of 2. We will therefore use the esti- 
mates in eg . (|33|) for the kaon and the eta, instead of dMi/drh = N^. Those are the 
states that will contribute more at low temperatures. For the rest, we will assume 
the uncertainty in dMi/drh = Ni to be from Ni/2 to 2Ni. 

Thus, in Fig. 4 we show the results when the massive states are taken into account. 
We have considered in eq. ( pTD all particles containing u and d quarks up to 1300 MeV 
and we have taken the central values of all the other parameters. The dominant 
contributions are, of course, those of the kaons, the eta, the rho and the omega. The 
shaded areas cover the uncertainty in dMi/dm that we have just described. Obviously, 
the net effect is biggest for the standard scenario, since the critical temperature is 
higher, where Tc is decreased down to 190 to 200 MeV. This result, although it has 
been obtained within the generalized formalism, reproduces very nicely the standard 
ChPT estimate given in 

Indeed, the r dependence is given in Fig. 5 where we plot the evolution of the 
chiral condensate both with T and r, for the central values of all the parameters, but 
also including the contributions from massive states. Note that, for the extreme case 
when r = 8, the decrease is of the order of 5 MeV, down to around 125 MeV. 
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Figure 5: The evolution of the chiral condensate when we include corrections in the 
pion gas from heavier states, as a function of the temperature and the quark mass 
ratio r. (Using the central values of all parameters and estimates) 

6 Conclusions 

In this work we have studied the generalized scenario of chiral symmetry breaking, 
either with a large or a small T = condensate. For that purpose we have described 
a pion gas by means of the virial expansion, whose coefficients have been calculated 
using the amplitudes obtained within O(p^) Generalized Chiral Perturbation Theory. 

We have also added a crude estimate of contributions from particles heavier than 
the pion, in a free gas approximation, which can be justified at low temperatures. The 
effect of these particles is always to decrease the temperature of chiral restoration. 
Their net effect is to lower Tc by 10 to 20 MeV in the standard scenario, and by 
around 5 MeV when the T = chiral condensate is smallest. 

From our results, it seems that the chiral phase transition in a pure pionic gas may 
occur at energies as low as 125 MeV in the lowest possible T = condensate scenario. 
The main source of uncertainty is the fact that within the Generalized approach many 
parameters still remain undetermined. In the worst case, which again corresponds 
to the lowest condensate and lowest Tc, it can be estimated at about 20%. For the 
standard case of a large condensate, we recover previous estimates of Tc — 190 MeV. 

In conclusion, we have found that the value of ~ 190 MeV for the critical tem- 
perature obtained from standard Chiral Perturbation theory can be seen as an upper 
bound if we were to include 0{m'^) corrections in the mass terms, in addition to the 
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standard condensate contribution. The effects of these corrections always lower the 
critical temperature, which, all together, could be as low as 125 MeV with a 20% 
uncertainty for the lowest condensate scenario. 
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